Coupled-Channel and Screening Effects in Charmonium Spectrum 
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Using the same quenched limit as input, we compare the charmonium spectra predicted by the 
coupled- channel model and the screened potential model in the mass region below 4 GeV, in which 
the contributions from decay channels involving P-wave (as well as even higher excited) D mesons 
can be neglected. We find that the two models have similar global features in describing charmonium 
spectrum since they embody the same effect of the vacuum polarization of dynamical light quark 
pairs. This agreement is encouraging for using the screened potential to further describe the mass 
spectrum of higher charmonium in the region above 4 GeV. Applying both models will be helpful 
to clarify the nature of the newly discovered "X, Y, Z" mesons in B-factories and other experiments. 
In particular, we show the S-wave decay coupling effect on lowering the Xci(2P) mass towards the 
DD* threshold, in support of the assignment of the X(3872) as a Xci(2P)-dominated charmonium 
state. 

PACS numbers: 12.39.Jh, 13.20.Gd, 14.40.Gx 



I. INTRODUCTION 



Studies on heavy qtiarkonium spectroscopy have been 
stimulated greatly in recent years by the discovery 
of many hidden charm states, the so-called "X, F, Z" 
mesons [l| in B-factories and other experiments. The 
QCD-inspired interquark potential models, such as the 
Cornell model 0, are successful in predictions of char- 
monium and bottomonium spectra below the open fla- 
vor thresholds. However, the existence of open charm 
thresholds can change the charmonium spectrum signif- 
icantly through the virtual charm meson loops. These 
coupled-channel effects were considered also in the Cor- 
nell model Pi , and techniques were further developed by 
the unitaritized quark model Q based on the ^Pq quark 
pair creation model Along this line and with the 
updated parameters, the coupled-channel effects in the 
charmonium spectrum have been further studied during 
these years These studies provide important in- 

formation on the identifications of the "X, Y, Z" mesons. 

The quark potential model is subject to modification 
due to quantum fluctuation, i.e., the creation of light 
quark pairs, which may be compensated by the virtual 
hadron loops in the coupled-channel model. For this rea- 
son, the quark potential model 0, which incorporates a 
Coulomb term at short distances and a linear confining 
term at large distances, may be called as the quenched 
potential model. On the other hand, the vacuum po- 
larization effect of the dynamical fermions may soften 
the linear potential at long distances Q, and cause the 
screening effect, which may be discussed phenomenologi- 
cally as the screened potential model 9., J^, . Such 
screening or string breaking effects have been demon- 
strated, although indirectly, by the simulations of un- 
quenched lattice QCD JJJ. This effect is also implied by 
calculations within some holographic QCD models [13| . 
The screened potential model has been used to reexamine 



the charmonium spectrum recently, and it is found 
that the masses of higher charmonium states are lowered, 
compared to the quenched potential model (Tsj , and the 
mass suppression tends to be strengthened from lower 
levels to higher ones. Such tendency can also be found 
in the calculations of the coupled-channels model, such 
as those in Rcf. 6). It is not very surprising since, in the 
parton-hadron duality picture, the two models embody 
the same effects of light quark pairs. 



Comparing to the screened potential model, the 
coupled-channel model is more difficult to handle in prac- 
tice, especially in the case when the P-wave ( and higher 
excited) D and Dg mesons as the intermediate states 
are involved. However, the latter can describe the near- 
threshold effect 0, UBl J which has been ignored in the 
former. It is then interesting to compare the two models 
using the same quenched limit as input in the domain 
of charmonium spectrum. The comparison has twofold 
meaning: the coupled-channel model can be helpful to 
establish the form of the screened potential and to deter- 
mine the screening parameter in the mass region below 
4 GeV, where the P-wave D mesons are expected to be 
decoupled; whereas the screened potential model can be 
helpful to normalize the global features of the coupled 
channel model in the mass region above 4 GeV. 



In this paper, we will compare these two models in 
the mass region below 4 GeV for charmonium spectrum 
with the same quenched limit. We will introduce the two 
models in turn. And we will compare results of these two 
models numerically, and finally a summary will be given. 
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II. QUENCHED AND SCREENED POTENTIAL 
MODELS 

We choose the CorneU model J^l as the quenched hmit, 
in which the potential has the well known form: 



V{r) 



h Ar 

3 r 



(1) 



where the first term denotes the color Coulomb force in 
the one-gluon exchange approximation due to asymptotic 
freedom of QCD at short distances, and the second term 
denotes the linear confining potential, which is consis- 
tent with both the rotating string picture [l7j and the 
quenched lattice calculations (see [1^ for a review and 
references) . The constant C in ([T]) is the renormolization 
term. 

In principle, one can relate the parameter etc to the 
running coupling constant as{mcv) in QCD, where « 
0.3 is the charm quark velocity squared in the char- 
monium rest frame, and relate A to the string tension 
T = l/(27ra') - 0.18 GeV^, where a' denotes the Regge 
slope in the rotating string picture Thus, we choose 



0.55, A = 0.175 GeV-^ 



(2) 



in the following analysis. 

To restore the hyperfine and fine structures of the 
charmonium spectrum, one needs to introduce the spin- 
dependent potential Vsdi which is relativistically sup- 
pressed compared to ^ in ([T]). Assuming the Lorentz 
structure of the linear confining force is of scalar type, 
the spin-dependent potential can be derived from ([1]) by 
the standard Breit-Fermi expansion to order jt? , and 
has the form [15j 
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where S = Sc + Se is the total spin, rric the charm 
quark mass, and the smeared delta function is taken 
to be S{r) = (cr/^F)^ e"'^'''' with a = 1.45 GeV 
Then, the Hamiltonian for the quenched potential model 
is given by 
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FIG. 1: Comparison of the Cornell potential V{r) (dashed 
line) and the screened potential T4cr (solid line) with r = 0.1-4 
fm. The constant terms are neglected and asymptotic limit 
of Vscr is shown in the dotted line. 



to substitute V{r). The first term on the right-hand side 
of ([5]) is taken to be the same as that in ^ due to its short 
distance nature. The screening parameter fj, sets the scale 
of distance at which the string breaks. We choose 



fi = 0.075 GeV, 



(6) 



of which the inverse is about two times of the radius of D 
meson. Needless to say, the right hand side of ([3]) should 
also be modified accordingly. 

At short distances, the screened potential Vscr in © 
has the limit consistent with that of V{r) in ([T|). How- 
ever, the difference between the potentials increases when 
r increases, and finally Vscr goes to a constant A/^ -I- C 
when r —^ oo, as shown in Fig. [TJ As a result, the char- 
monium spectrum in the screened potential model will be 
suppressed compared to that in the quenched potential 
model, especially for the higher exited states. Further- 
more, in the screened charmonium spectrum, there will 
be a saturation energy of about 5-6 GeV, at which the cc 
quark pair can not be bound together at all. 

Let's stress again that the parameters in ^ and ^ 
may not be the same as those chosen in [l^l, since our 
purpose here is to compare the spectra in the two mod- 
els rather than to get a good fit of the spectrum to the 
experimental data. 



Ho = 2mc 



+ Vir) + Vsd{r), 



(4) 



where the kinematic energy term have been included ex- 
plicitly. 

As mentioned above, the linear confining potential will 
be softened by the vacuum polarization induced by the 
dynamical quark pair creation. Such unquenched effect 
can be roughly accounted for by modifying the l ong dis- 
tant behavior of V{r) in Q. Following Refs. [l3,lll|, we 
use the screened potential 



Vscr{r) 
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III. 



THE ^Po-BASED NON-RELATIVISTIC 
COUPLED CHANNEL MODEL 



The light quark pair creation from vacuum is assumed 
to share the same quantum number, 0^^, as the vacuum 
in the '^Pq model [4]. In the non-relativistic limit, the 
■^Pq model can be represented by the Hamiltonian [19l | 



QPC 



(7) 



where niq is the mass of the produced quark, and 7 
is a constant reflecting quark pair creation strength. 
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which can be determined phenomenologically by the OZI- 
allowed decay widths of charmonia There are no 

fundamental reasons to enhance the heavier quark cre- 
ation relative to lighter ones, which is implied by the 
factor of 2mq in ([7]). Therefore, following Ref. Q, we use 
the effective strength 7s — ^7 for the strange quark, 
where ^ is the ratio of constituent masses of up and 
down quarks(mo = m„ = nid) and strange quark(r7is). 
Following Ref. Q, we choose 7 = 0.322. 

Light quark pair creation can result in mixing between 
the bare charmonium state ijjQ{cc) and the open charmed 
meson pair B{cq) and C(qc). Thus, neglecting the mixing 
among the bare states the physical state 5" can be 
represented by 



I = ao I ^'0) + 



BC 



J dv CBc{i^) I BC;v) 



(8) 



where u denotes the variable of three-momenta of the BC 
system. The coefficients ao and cbc{i^) are understood 
to be subject to the normalization of the corresponding 
wave functions. 

The Hamiltonian of this coupled-channel system can 
be formally written as 



BC 



QPC- 



(9) 



Here, Ho has been introduced in ^ but can be different 
from the one in the quenched potential model by different 
renormalization constant C. The Hamiltonians Ho and 
Hbc can only act on the states -00 and BC, respectively, 
and give the bare spectra of them: 

Ho I V^o) = Mo I Vo), (10) 
Hbc I BC) = EBcity) I BC;i^), (11) 

where Mq is the so-called bare mass of the bare charmo- 
nium state and 



Ebc{Pb,Pc) = +P^ + ^M^ + P^, (12) 

provided that the inner interactions between B and C 
can be neglected. On the other hand, the Hqpc in ([9]), 
which is defined in ([7|), can acts only between | ^0) and 
\BC-v). 

The Hamiltonian H in ([9]) defines the physical mass M 
of the state 5' as 



I = M I 



(13) 



and the mass M can be obtained by solving the multi- 
channel Shrodinger equation Substituting ((5|)-pT|) 
into (|13p. one can get the integral equation 



M + n(Af ) - Mo = 0, (14) 
where the self-energy function H{M) is given by 
\{BC-v\Hqpc\^o)? 

BC •' 



n(M) ^Y. h"" 



Ebc{v) - M + it 



(15) 



When Re[M] > Ma + Mb, it is obvious that the function 
n(Af ) will not be real, and the imaginary part is propor- 
tional to the decay width of ^' ^ BC. Therefore, one can 
solve Eq. (fH|) in the complex plane and get the pole mass 
Re[M] [6,]. However, we will define the coupled-channel 
mass Mcou as 



M,o„ + Re[H(Af,„„)] - Mo = 0, 



(16) 



which is also called the Brcit-Wigner mass by the authors 
of Ref. ;6]. 

It is worth emphasizing here the difference between 
the definitions of the bare mass in and of ours. The 
authors of Ref. [1] once subtract the dispersion integral 
in (fT5|) at Ad^,, and absorb the term H(M^) into the bare 
mass definition: 



=Mo-H(M^), 



(17) 



where H(M^) is real and positive, which can be seen di- 
rectly from the definition of the function H(M) in (|15p . 
If the subtracted constant n„(M^) are the same for all 
the charmonium states n, the bare mass Afp in (|17p is 
just a rescaling one of Mo and the renormalized mass 
shift would not changed. In Ref. the matrix element 
square \{BC \ Hqpc \ V'o)!^ in CSl) is simply parame- 
terized by using a exponential form factor, and the node 
structure in the wave function of higher excited state 
is absolutely neglected. Thus, the subtracted constants 
H„(Af^) for these excited states might be overestimated, 
and as a result, the renormalized mass shifts in Ref. Q 
are commonly smaller than those in our model, as one 
can seen in the following section. 

On the other hand, the wave functions of "00: B and 
C are needed to determine the matrix element {BC \ 
Hqpc I "00 ) in These wave functions are usually 

chosen as the harmonic oscillator ones [1, 0, • How- 
ever, we determine them by solving the quenched mass 
equation Eq. (fTO|) with Eq.Q. But for simplicity, we will 
neglect the corrections due to the spin-dependent poten- 
tial Vsd to these wave function. 



IV. NUMERICAL RESULTS AND DISCUSSION 

As mentioned in Sec. I, our aim is to compare the 
charmonium spectra or the renormalized mass shifts 



AMcou = Mc. 



liquet ^^Iscr — 



A'Ve, (18) 



where the subscripts que, cou, and scr denote the 
results obtained from the quenched potential model, 
the coupled-channels model, and the screened potential 
model, respectively. The quenched mass Mq^e can be re- 
lated to the bare mass Mo in the coupled-channels model 
by the relation: 



M„. 



Co 



Mo - a, 



(19) 



4 



TABLE I; Chamonium spectra and mass shifts in different models in units of MeV. Here, the subscripts que, cou, and scr 
denote the results obtained from the quenched potential model, the coupled channel model, and the screened potential model, 
respectively. The mass shifts AMcou and AMscr are listed in the 5th and 6th columns, respectively. The results of Ref. [9] are 
also listed. The bare mass in the 7th column is copied from Ref. p^ . All the quantities listed here should be understood as 
the renormalized ones. 



states 


Our results 


Results of Ref. [6] 




Mean 


Mscr 


AMcou 


AMscr 


Mo 


M' 

^'-^cou 




— -. 1 ^ 

1^5*0 


2980 


2980 


2980.0 








2982 


2982 





TTT^ 


3112 


3100 


3105 


-12 


-7 


3090 


3090 





^ 1 


3583 


3531 


3539 


-52 


-44 


3516 


3514 


-2 


^ '-i 

I'^Po 


3476 


3441 


3448 


-35 


-28 


3424 


3415 


-9 


I'' Pi 


3568 


3520 


3526 


-48 


-42 


3505 


3489 


-16 




3628 


3565 


3577 


-63 


-51 


3556 


3550 


-6 


2' So 


3697 


3635 


3626 


-62 


-71 


3630 


3620 


-10 




3754 


3674 


3674 


-80 


-80 


3672 


3663 


-9 




3895 


3818 


3805 


-77 


-90 


3799 








3878 


3794 


3790 


-84 


-88 


3785 


3745 


-40 


i^D2 


3896 


3818 


3805 


-78 


-91 


3800 








3903 


3823 


3812 


-80 


-91 


3806 






2' Pi 


4042 


3961 


3909 


-81 


-133 


3934 


3929 


-5 


2'Po 


3948 


3915 


3839 


-33 


-109 


3852 


3782 


-70 


2^ Pi 


4030 


3875 


3900 


-155 


-130 


3925 


3859 


-66 


2^ Pa 


4085 


3966 


3941 


-119 


-144 


3972 


3917 


-55 



where Cque and Ccou denote the renormalization con- 
stants in the quenched potential model and the coupled- 
channels model, respectively. Thus from Eq. (|16p . the 
mass shift AMcou can be given by 

AMcou = -Re[n(Meo„)] - Cgue + C^ou- (20) 

To improve the reliability of the calculation of the 
coupled-channel model, we restrict the mass region of 
charmonium spectrum to be below 4 GeV, in which the 
decays to P-wave ( and higher excited ) D and Dg mesons 
are kinematically forbidden, and only the S-wave D and 
Ds mesons are involved. So in the following we use the 
S-wave mesons, of which the masses are well de- 
termined and the widths can be neglected [23|, as the 
intermediate states only. 

In the mass region below 4 GeV, the charmonium spec- 
trum consists of the IS, 2S, IP, 2P and ID levels. For con- 
venience, we choose the renormalization condition such 
that the predicted tJc masses in the three models are fixed 
to be the measured value [13], i.e., 2980 MeV. Together 
with the inputs of the quark masses 

rUc = 1.7 GeV, to, = 0.33 GeV, m, = 0.5 GeV, (21) 

the renormalization condition determines the constant 
terms in dT]): 

Cgue = -419 MeV, Ccou = -272 MeV, C^cr = -403 MeV. 

(22) 

We list the numerical results of the mass spectra and 
the mass shifts in Tab. HI The subscripts que, cou, and 
scr denote the results obtained from the quenched poten- 
tial model, the coupled channel model, and the screened 



potential model, respectively. And the mass shifts have 
been defined in ((THl). In Tab.[Tl both AMcou and AMscr 
have minus sign, and on the whole, they are consistent 
with each other. This can also be seen from Fig. 2. 

In addition, for the IS, 2S, IP and ID levels, which 
either couple the D meson pair in P-wave or lie far away 
from the threshold of the D meson pair, the mass shifts 
AMcou s are comparable to each other, which is consis- 
tent with the first Hadron Loop Theorem derived by the 
authors of Ref. [5] in the approximation of equal masses 
for charmed mesons. 

However, if the cc pair couples to D meson pair in S- 
wave and the mass M calculated in the coupled channel 
model (namely Mcou) in (|15p is close to the threshold 
Mb + Mc, the self-energy 11, then the mass shift AMcou, 
will strongly depend on M{Mcou)- This is well known as 
the S-wave threshold effect [l^, and as a result, the loop 
theorem will be violated. 

This is just the case for the 2P charmonium states, 
since for some of them the coupled-channel masses Mcou 
(see Tab.|l]) are close to the thresholds of the S-wave chan- 
nels DD* + c.c. and D*D*^. More in details, for the 2^Po 
state (Xco)' coupled-channel mass Mcouix'co) — 3915 
MeV, which is fairy far from the threshold of the S-wave 
channels DD and D*D* . Consequently, the mass shift is 
as small as 33 MeV. As a second example, the coupled- 



The contributions to the mass shifts from the Dg* Dg* chan- 
nels are very small due to the small strange quark pair creation 
strength 73. 
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FIG. 2; Mass shifts AMco« and AMs, 
mass AIq. 



varying with bare 



FIG. 3: The physical mass Mcou dependence of the unrenor- 
mahzed mass shifts — Re[n(Mcou)] for the 2^ Pi state. 



channel masses of the 2^Pi [h'^) and 2^P2 (Xc2) states are 
roughly equal. However, their mass shifts induced by the 
D*D* , of which the threshold is closest to their coupled- 
channel masses, are different by a factor of 2 (see Tab. 
I in Ref. 0). As a result, the mass shift of 2^Pi state 
is smaller than that of 2^P2 state. Finally, the coupled- 
channel effect of (Xci) state should be most signifi- 
cant since the mass Mcou{Xci) = 3875 MeV is very close 
to the threshold of D°D*° /D+D*- +c.c.. This result can 
also give support to the x'd assignment of X{3872) [2l| . 

It needs emphasizing here that the closeness of 
Mcouix'ci) to the threshold of D°D*° is not very sen- 
sitive to the bare mass of x'd- This can be seen from 
Fig. 3, where the physical mass Mcou dependence of the 
unrenormalized mass shift —Ke[Il{Mcou)] for the 2'^Pi 
state is shown. The relation between the unrenormalized 
mass shift and the renormalized one AMcou is given in 
(pO)) . From Fig. 3 one can see the mass shift function is 
strongly dependent on the physical mass. As a result, 
the slope of the mass shift curve is very large near the 
threshold and "attract" the mass Mcouix'ci) towards the 
threshold. That means, if one changes the bare mass, 
say, by 50 MeV, the change of the mass Mcouix'ci) is 
only about 10-15 MeV as can be roughly found from the 
figure. This is just a realization of the S-wave thresh- 
old effect [IQ] in our coupled-channels model. More pre- 
cisely, the curve in Fig. 3 shows the cusps in the neutral 
and charged D*D channels numerically (see the second 
paper in [lg| for more discussions). Thus, the physical 
mass of x'ci is quite natural to be close to the threshold 
of D^D*^ and thus x'd "^^Y be a good candidate for the 
A:(3872) f2ll. 

We also list the results of Ref. @ in Tab. U for com- 
parison. As we have mentioned, the definition of the 
renormalized bare mass Mq ([T7]) is a little different from 



that of Mq. Moreover, the node effect in the wave- 
function overlap integral for a higher excited charmo- 
nium decaying into charmed mesons is not considered 
in Ref. Q. As a result, the renomalized mass shift 
AM^„„ = n(M^) - Ren(M) in ^ tends to be smaller 
than the one in (pO)) for excited states as has been an- 
alyzed in the last section. Furthermore, the authors of 
Ref. @ use the available results of the quenched poten- 
tial model [IBl, in which different parameters from ours 
are chosen, as their bare mass inputs (the 7th column in 
Tab. [l|. However, comparing the mass shifts in the 5th 
and 9th columns in Tab. HI one can find that they both 
have similar features as that of AMgcr- This indicates 
that the screened potential in ^ depicts the main fea- 
ture of the vacuum polarization effect of the dynamical 
quark pair creation, although it fails to describe some 
fine structures, such as those induced by the near S-wave 
threshold effects. 

To compare the global features of spectra for the 
coupled-channel model and the screened potential model 
more directly, we also illustrate the c.o.g. shifts AM^°^ 
and AM™^ in Fig. 4. Here, the c.o.g. is defined as 



]\^cog 



M,, 



3M-^ 



(23) 



where Mgig denotes the mass of the spin-singlet state, 
and M^°f the c.o.g. of the spin-triplet states. From 
Fig. 4, one can see both the mass shift trajectories ex- 
hibit good behaviors, and they are roughly consistent 
with each other as the above analysis. However, the in- 
crease of the mass shift AMj°? tends to be faster than 
that of AMcou iu the higher mass region. This seems to 
indicate that the potential in ([5]) somewhat overestimates 
the screening effect. But it is not the whole story since 
the P-wave ( and higher excited ) D(^s-^ mesons contribu- 
tions, which have been neglected, also tends to enhance 



6 



20- IS 



2P 



' cou 

|AM™3 
' scr 



ID 



2S 
* 



IP 



3200 3400 3600 

Bare Mass (MeV) 



FIG. 4: The trajectories of the c.o.g. shifts AM^°^ and AM^°^ 
varying with the bare mass Mq. 



in the mass region below 4 GeV, in which the contribu- 
tions from channels involving the P-wave ( and higher 
excited )D and Ds mesons can be neglected. We use 
the same quenched limit for the two models. And for 
the coupled-channel model, we use the wave functions 
obtained by the quenched potential model in the non- 
relativistic limit to determine the hadronic transition ma- 
trix elements (wavefunction overlap integrals) in psp . 

We find that the two models have similar global fea- 
tures in describing the charmonium spectrum as expected 
since they embody the same effect of the vacuum polar- 
ization of the dynamical quark pair creation. This agree- 
ment is encouraging for using the screened potential to 
further describe the mass spectrum of higher charmo- 
nium in the region above 4 GeV. Applications of both 
models will be helpful to clarify the nature of the newly 
discovered "AT, y, Z" mesons 1]. In particular, we show 
the S-wave coupling effect on lowering the Xci i'^P) mass 
towards the DD* threshold as support to assign the 
X(3872) as a Xci(2P)-dominated charmonium state. 



the mass shift AM^°^ in the same higher mass region. 

V. SUMMARY 

We compare the charmonium spectra predicted by the 
coupled-channel model and the screened potential model 
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